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TO ALL WHOM IT MAY CONCERN: 

BE IT KNOWN that we, Hiroo OHMORI , a subject of Japan 
and residing at Yokohama-shi, Kanagawa, Japan, Takahiro ASAI a subject of 
Japan and residing at Yokosuka-shi, Kanagawa, Japan and Tadashi 
MATSUMOTO, a subject of Japan and residing at Ebina-shi, Kanagawa, 
Japan have invented certain new and useful improvements in 

"ADAPTIVE EQUALIZATION METHOD AND ADAPTIVE 

EQUALIZER" 

and we do hereby declare that the following is a full, clear and exact 
description of the same; reference being had to the accompanying drawings 
and the numerals of reference marked thereon, which form a part of this 
specification. 
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ADAPTIVE EQUALIZATION METHOD AND ADAPTIVE EQUALIZER 

BACKGROUND OF THE INVENTION 

The present invention relates to an adaptive equalization method and 
5 an adaptive equalizer for use, for example, in mobile communications to 
reduce the influence of fading. 

In a mobile communications environment signals transmitted from 
mobile stations arrive at a base station from various directions under the 
influence of surrounding buildings, trees and so forth. The influence of 
1 0 intersymbol interference, which results from the arrival of paths of large delay 
time differences in such a multipath channel, poses a serious problem for 
mobile communications. A known technique for lessening the influence of a 
signal waveform distortion by the intersymbol interference is an equalization 
technique that estimates the characteristic of the channel and performs 
1 5 equalization of delayed paths based on the estimated characteristic. One of 
equalization algorithms used in the equalization technique is a 
turbo^c oding-technique-ag plied iterative equalization algorithm that achieves 
ayery excellent err or rate characteristic in low SNR, The turbo coding 
method implements, with practical coding-decoding complexity, an error rate 
20 characteristic that is almost a ceiling value obtainable with a channel coding 
method; only the presence of such an error rate is merely suggested in 
information theory. With the turbo coding method, the transmitting side 
transmits coded information sequence and its rearranged or interleaved 
version, and the receiving side repeatedly decodes them alternately with each 
25 other. Sin ce the process p f generating a coded sequence of convolution 
codes in the turbo coding is equivalent t o the process of generation of the 
^iiilersymboi^^ in multipath channel, the decoding process of the 
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turbo coding technique can be applied to the equaliz ation ofd elaved paths. 

Fig. 1 illustrates in block form an example of a transmitter-receiver for 
the iterative equalization disclosed in Daryl Reynolds and Xiaodong Wang, 
"Low Complexity Turbo-Equalization for Diversity Channels," 
5 http://ee.tamu.eduT revnolds/ . A transmitter 10 encodes an information bit 
sequence d(i) by an encoder 11, then rearranges the coded output sequence by 
an interleaver 12, and modulates the carrier by a modulator 13, thereafter 
transmitting the modulated output bit sequence b(i) over a channel 14. In a 
receiver 20 delayed paths are equalized by a SISO (Soft-Input Soft-Output) 
|* 10 equalizer 21 for the received signal r(t), and the log-likelihood ratio (LLR) Ai 
|3 between +1 and -1 for each code bit b(i) of the received signal r(t) is derived 

by the following equation (1). 
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11 KJ1 B Pr[b(i) = -l|r(t)] 



= wM^l + lQ Pr[b(i) = + 1] 
|J *p[r(t)|b(i) = -l] *Pr[b(i) = -l] 

b + Mb(i)] (i) 



In Eq. (1) the first term, Aa[b(i)], is extrinsic or external information 
15 derived by subtracting a priori information A, 2 [b(i)], by a subtracter 26, from 
the log-likelihood ratio A] [b(i)] provided from the SISO equalizer 2 1 . The 
second term, A, 2 [b(i)], is the a priori information that is a rearranged or 
interleaved version, by an interleaver 25, of X 2 [b(n)] obtained by subtracting 
Xi[b(n)], by a subtracter 24, from a log-likelihood ratio A 2 [b(n)] provided 
20 from a SISO decoder 23 in the previous iteration. The initial value of 

A, 2 [b(n)] is zero. Th^priori information A, 2 [b(i)] is used as an a priori LLR 
value of a coded bit string {b(i)}. The extrinsic information is restored by a 
de-interleaver 22 to Ai[b(i)], which is input as a posteriori information to the 
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SISO decoder 23. The SISO decoder 23 uses th^j^osteriori information to 
derive a a posteriori LLR value A 2 of the coded bit string (b(i)} by the 
following equation (2). 

A 2 [b(n)]=^[b(n)J+^[b(n)] (2) 
5 In the second and subsequent iterations, Xi[b(n)] is subtracted by the 

subtracter 24 from the a posteriori LLR value A 2 derived by the SISO decoder 
23 to obtain extrinsic information X 2 [b(n)], which is provided as a priori 
information A, 2 [b(i)]via an interleaver 25 to the SISO equalizer 2 1 . Such 
iterations of equalization and decodi ng update the LLR valu e, providingj_ 
|* 10 more probable LLR value. That is, the first log-likelihood ratio Ai is 
j«j calculated by the first adaptive equalization of the received digital signal, and 

*2 the first log-likelihood ratio Ai is decoded to obtain the second log-likelihood 

in 

) Q ratio A 2 , which is used to perform the second adaptive equalization of the 

2 received digital signal to obtain the first log-likelihood ratio Aj, which is then 

I- 

f U 15 decoded, and the decoded value is output. 

I* 

(0 In the iterative equalization the SISO equalizer and the SISO decoder 

\1 use a MAP (Maximum A-posteriori Probability) algorithm, Log-MAP 

algorithm based on the MAP algorithm, or Max-Log-MAP algorithm that has 
reduced the computational complexity of equalization and decoding through 
20 use of approximate calculation. Still another algorithm is S OVA 

(Soft-Output Viterbi Algorithm) that provides bit reliability information while 



performing opt imum dec od ing o n a sequence-wise basis. In the case of 
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using any of these algorithms for equalization by the equalizer, however, the 
complexity of the algorithm exponentially increases with respect to the 
25 maximum number of delayed symbols to be considered in the equalizer, 

resulting in very high computational complexity. For example, when BPSK 
is used as the modulation system, assuming that the storage capacity of one 



receiving channel is 10 symbols (the maximum number of delayed symbols to 
be considered in the equalizer being 10 symbols), the number of states 
becomes as large as 2 10 =1024, exceeding practical computational complexity. 
The computational complexity could be reduced by using a linear equalizer 
5 for the iterative equalization. 

Fig. 2 illustrates a typical example of the linear equalizer. In the 
illustrated equalizer sampled received signals are applied to a 
^ascade -conne ction of delay elements (of x-sec delay) 3 1 to obtain signals 
delayed by x-sec intervals in a sequen tial order, which are weighted by a 
1 0 variable weighting circuit 32 having set therein a coefficient by a coefficient 
adaptation circuit 34 and combined by an add er 33 to implement equalization. 
As the delay time x, the sampling perjod Xor T/2J susually chosen—- 
Considering me eq^trafeEraT one filter, the linear equalizer achieves 
equalization by implementing a filter characteristic convenient for identifying 
1 5 the received signal. This filter basically implements an inverse characteristic 
of a channel distortion. 

Wheiisuch a lmear^ializer is used for the abo ve -described itera tive 
e qualiza tion, it is necessary to derive a soft decision value as the equalizer 
output. In the afore-mentioned literature by Daryl Reynolds and Xiaodong 
20 Wang, there is set forth a scheme in which, by deriving the soft decision value 
through the use of the linear equalizer, the iterative equalization can be 
achieved with less computational complexity than in the case of using the 
conventional algorithm such as MAP, SOVA or the like. In the literature a 
description is given of the application of the linear equalizer assuming the 
25 following discrete-time model as a channel model. A received signal sample 
r(n) is expressed as follows: 
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Kn) = Ib(n - j)h(n; j) + ov(n) , n=0, 1, N-l (3) 
j=o 

v(n) =[v 0 (n) v,(n) . . .v M .,(n)] T (4) 
h(nu) =[ho(n;j) h,(n;j) . . .h M -i(n:j)] T (5) 
In the above, M is the number of channel outputs and 
r(n)=[ro(n)r,(n)...r M -i(n)] T is a vector representing a received sample value in 
each channel at time n, where T indicates a transposed matrix. For example, 
in the case of using a plurality of antennas, M corresponds to the number of 
antennas. B(n) represents a coded bit, and v(n) in Eq. (4) is a complex 
vector of an average 0, representing noise. In Eq. (5), h(n j) is a vector 
representing the tap weighting coefficient of the channel, and J is the channel 
ISI length. Now, a received signal vector R(n), a noise vector V(n), a coded 
bit vector B(n) and a channel matrix H m (n) are defined as follows: 
~r(n + J-l)" 

R(«)= . : _ (6) 

-MJxl 



V(n) = 



r(n) 
v(n + J-l) 

v(n) 



(7) 
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15 B(n)=[b(n+(J-l))b(n+(J-2))...b(n+l)b(n)b(n-l)...b(n-(J-l))] T 
H(n) = [h(n;0)-h(n;J-l)] 

"H(n) 0 0 

0 H(n) i 

: *•. 0 

0 •• 0 HCn) 



(8) 
(9) 
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h(n;0) h(n;l) ... h(n;J-l) 0 
0 h(n;0) h(n;l) ••• h(n;J-l) 



0 
0 



h(n;0) h(n;l) 
0 h(n;0) 



h(n;l) 



0 

0 0 

h(n;J-l) 0 

h(n;J-l) 
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The vectors R(n) and V(n) given by Eqs. (6) and (7) are column vectors 
of a matrix with MJ rows and one column, that is, MJ-dimensional column 
vectors. B(n) is a (2J-1 )-dimensional column vector, H(n) is a matrix with 
M rows and J columns, and H m (n) is a matrix with MJ rows and (2J-1) 
columns. From the above, Eq. (3) becomes as follows: 
R(n) -HnOi^n) +oV(n) (1 1) 

Next, the likelihood of the code bit string is defined by the following 
equation using the a priori LLR value A, 2 [b(k)] of the code bit string b(i). 
/ b»(k)=tanh[A 2 [b(k)]/2], n-(J-l) <k<n+(J-l) (12) 
Eq. (12) has the property that it rapidly approaches 1 or -1 upon 
deviation of a variable x of a hyperbolic function tanhx from 0, and the 
likelihood b' (k) represents an estimated value of b(k). Further, that 
estimated value vector of the code bit in the received signal vector R(n) which 
affects, as intersymbol interference, code bit b(n) at time n is given by the 
following equation. 



B , (n)=[bXn+(J-l))b , (n+(J-2))...b'(n+l)0b'(n-l)...b'(n-(J-l))] 1 



(13) 



In Eq. (13), since F(n) is an estimated value vector that is used to subtract 
fromjthereceived sig naFvector Rfc Qthe intersytnborinterference component 
fo^^cocLe-bit b(n) jtthrietuthe element corresponding to b(n) of the vector 
-K(n)js set to zero. A difference vector Rg(n) by subtracting an intersymbol 
interference component replica from the received signal vector R(n) is 
defined by the following equation using the estimated intersymbol 
interference code vector B'(n). 



Rc(n) = R(n)- H m (n)B r (n) = H ra (n)[B(n)-B , (n)]+av(n) (14) 
The difference vector Rc(n) is subjected to filtering expressed by the 
following equation through the use of a linear filter of tap coefficient G(n). 

Z(n) = G(n) H R c (n) (15) 
where H is a conjugate transposed matrix. The tap coefficient G(n) is so 
^choseiias jo minim izeti ie mean square error between the code bit B(n) and 

the filter outputZ(nVa^expressed h y th fi follo w in g equation 

G(n)=arg min Efllb^-GfcfR^n)!! 2 } 

=arg min G(n) H E{R c (n)R <; (n) H }G(n)-2 G(n) H E{b(n)R c (n)} (16) 
£ 10 where arg min E{} represents that one of all J-dimensional vectors G whose 
|j content of { } is minimum. 

1^ Expansion of first and second terms in Eq. (16) gives 

g E{R.(n)R c (n) H }=H m (n)A(n)H m (n) H +a 2 I (17) 

| 4 E{b(n)R c (n)}=H m (n)e J (18) 

j U 15 where I is a unit matrix and 

«j A(n)=Cov{B(n)-B'(n)} 

U ===diag[l-b ,2 (n-KJ-l)),... J l-b'Vl),l,l-b ,2 (n-l), 

...,l-b' 2 (n-(J-l))] (19) 
In the above, Cov represents a covariance matrix and diag represents a 
20 diagonal matrix. e ; is a vector of a length 2J-1 , a J-th element being 1 and 
other elements being Os. By solving the normal equation with the gradient 
vector on the right-hand side of Eq. (16) to 0 and substituting Eqs. (17) and 
(18), the optimum tap coefficient G(n) is given as follows: 

^Hft*^"fri^>i (20) , 

25 Substitution of the thus derived tap coefficient G(n) into Eq. (1 5) gives 

Z(n)=e, T H m (n) H p m ^ (21) 
In this case, since it is necessary that H m (n) be a full column rank, the antenna 



number M needs to be at least 2 or more. 

Assuming that the filter output can be approximated by a Gaussian 
distribution, Z(n) by Eq. (21) can be described as follows (see V. Poor and S. 
Verdu, "Probability of Error in MMSE Multiuser Detection" IEEE Trans. 
Information Theory., vol. IT-43, No. 3, pp. 858-871, May 1997). 

Z(n)=u(n)b(n)+n(n) (22) 
where u(n) is an equivalent amplitude of the output signal and r)(n) represents 
a Gaussian distribution of an average 0 and the variance v^n). Therefore, 
u(n) and v(n) can be expressed as follows: 

Kn)=E{Z(n)b(n)} 

=[H m (n) H [H m (n)A(n)H ro (n) H +o 2 I]- , H m (n)3 J>J (23) 
where indicates the element at the intersection of the J-th row and the J-th 
column of the matrix. 

v 2 =var {Z(n)}=u(n>n 2 (n) (24) 
From the above, the extrinsic information derived by the linear equalizer can 
be deduced from the following equation. 

X,[b(»)] =1 og «^ = 4^)} 
p(Z(n)|b(n) = -l) l-u(n) 

To obtain the optimum tap coefficient G(n) with the above-described 
method, however, it is necessary that the inverse matrix calculation expressed 
by the following equation be conducted in Eq. (20) for each point in time; 
hence, much time is required for calculation. 

<D(n) = [H m (n)A(n)Hn,(n) H + c 2 If l (26) 
With the conventional method for iterative equalization by the linear 
equalizer, the inverse matrix calculation needs to be conducted for each point 
in time to update the tap coefficient-this gives rise to the problem of high 
computational complexity. 
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SUMMARY OF THE INVENTION 

It is therefore an object of the present invention to provide an adaptive 
equalization method and apparatus that permit reduction of the computational 
complexity for deriving the tap coefficient of the linear equalizer for iterative 
equalization and hence allow shortening the time for computation. 

According to a first aspect of the present invention, aregdyelsjgnaj 
s^pjejaluTandT^wn signal areusTd to estimate the impulse response~of 



the c hanneyand a received signal sample value sequence is used to adaptively" 
calculatejhe tap coeffici ent of a linear filter. Th g tap coefficient is us ed to 
P 1 0 ^. perfam linear filterin g^fflie^cei yed signal sa ngle^ajue^equence, and the 
j J result of lmear filtering and Mises timated value of the imp ulse response are 

! " " " ' ~ — 1 

S ^ used toje ^yeasoft decisionvalue. 

jj[ With this construction, in the case of using the linear equalizer for 

* iterative equalization, t he received signal sam ple value sequence is input 

j U 15 thereto to adapti vely derive the tap coefficient oTthelmear hlter-thiTpgrmits 
J J red uction of compu tational compl exity as compare^wjtlrthe prior art method - " 

I* involving the mverse tmatrix calculation. 

According to a second aspect of the present invention, in second and 
subsequent processing by an adaptive iterative equalizer for the iteration of 
20 equalization and decoding, an est imated v alue vector of the code bit based on 
the likelihood of a soft decision value from a decoder is subjected to linear 
filtering with an estimated impulse resp onse value H m (n) to generate a replica 

*^ — — - ■ ■ . ■ " ■ . — 

ofthejwdyedsigna^ The re plica is subtracted from the received sig nal to 
cr^te^if&rencesignal Rc(n) having removed t herefrom the intersy mbol 
25 inter ference componen t, then t he tap coefficient of the linear fil ter is * ' * 
calculated, for example, using the Matrix Inversion Lemma, from the / ' 

estimated jmpulseresponse value vector H L (n) approximated with the 
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intersymb d interfer ence component regarded as zercQhen,the difference, 
then the difference signal Rc(n) is subjected to linear filteTmgusing the tap 
coefficient, and the result of linear filtering and the estimated impulse 
response value^vecto rare used to derive the soft decision value. 
5 With this construction, the tap coefficient of the linear filter in second 

and subsequent iterative equalization by the linear equalizer can be calculated, 
for example, using the Matrix Inversion Lemma, by which it is possible to 
achieve adaptive equalization with less computational complexity than in the 
case of the conventional that involves the inverse matrix calculation of Eq. 
j: 10 (20). 

li According to a third aspect of the present invention, a received signal 

sample value sequ ence and a known signal are used to es timate the im pulse i 
[J J^sp^ise^fth^diaimel, then a threshold value is set accordingly, and 
. ! . jcomp^nentsofVawes^iaall^^ value are removed from the 

asi * "> ^_ 

| U 15 estimated impulserg^ns e-s^alijevector, that is, the dimension of the 
jO estimatedj ni pulsej^spo me ^lue vector is reduced. Further, components 

I * corresponding to those o fthe values smaller than the threshold value and 

hence removed from the estimated impulse response value vector are removed 
from the tap coefficient ofa linear filter calculated from the received signal 
. 20 samplevalue s equence and the re ceived signal vector to thereby reduce their 
dimension as well. The soft decision value is derived using the estimated 



impulse res ponjejyahje_ygctor, th e tap coeffic ient ofthe linear filter and the 
r eceived signal vector of mej husreduced dimensions. 

The above scheme permits reduction of the computational complexity 
25 for deriving the soft decision value when the linear equalizer is used for 
iterative equalization. 

According to a fourth aspect of the present invention, the likelihood of 
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the soft dec isio n value from a decoder is calculated in seco nd and subsequent 
processing by an adaptive equalizer for iteration of equalization and decoding, 
then the tap coe fficient of the linear filter is calculated, for instance, using the 
, Matrix Liv ersion Lemma from the estimated impulse response value vector 
reduced in dimension in the first round of iteration and stored in a memory, 
then the tap coefficient is used to perform linear filtering of a difference signal 
having subtracted therefrom an intersymbol interference component, and the 
result of linear filtering and the estimated impulse response value vector are 
used to derive the soft decision value. 



[jj 1 0 This scheme allows reduction of the computational complexity for 

|= deriving the soft decision value in the second and subsequent iterative 

[jj equalization by the linear equalizer as well. 



BRIEF DESCRIPTION OF THE DRAWINGS 



[ U 15 Fig. 1 is a block diagram of a transmitter-receiver of the iterative 

10 equalization scheme; 

(3 

Fig. 2 is a block diagram of a linear adaptive equalizer; 
Fig. 3 is a block diagram illustrating the functional configuration for 
first equalization processing in a first embodiment of the adaptive equalizer 
20 according to the present invention; 

Fig. 4 is a block diagram of the functional configuration in second and 
subsequent rounds of equalization; 

Fig. 5 is a graph for explaining the effect of the present invention; 
Fig. 6 is a diagram depicting an example of a multipath channel 
25 environment; 

Fig. 7 is a block diagram illustrating the functional configuration for 
first equalization processing in a second embodiment of the adaptive equalizer 
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according to the present invention; 

Fig. 8 is a block diagram of the functional configuration in second and 
subsequent rounds of equalization; 

Fig. 9 is a graph showing an example of an impulse response 
contained in the received signal; 

Fig. 10 is a graph showing an example of selecting etfective paths 
from the impulse response; 

Fig. 1 1 is a graph for explaining the effect of the present invention; 

and 

Fig- 12 is a graph for explaining a channel model used for computer 
simulations. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
EMBODIMENT 1 

In the method set forth in the afore-mentioned literature that employs 
a linear equalizer for iterative equalization, the likelihood b'(k) of the code bit 
is calculated by Eq. (12) using the a priori LLR value X 2 [b(k)] of the code bit 
string. I n the first equalization processing, since the a priori LLR value 
hlbQQ] in Eq. (12) is zero in the absence of the a priori information, th e 
likelihood b'(k) is all zeros as shown below. 

b'(k)=[00...0...] T 
Therefore, the difference vector R«(n) by Eq. (12) becomes as follows: 

Rc(n) =R(n) -H m (n)B'(n) =R(n) (27) 
and A(n) becomes a (2J-l)x(2J-l) whose diagonal components are all ones 
as shown below. 

A(n) = Cov{B(n) -B'(n)}=diag{l,l,...,l 1,1} (28) 

Accordingly, the optimum tap coefficient G(n) of the linear filter by 
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Eq. (20) can be derived by the following calculation. 

G(n)= [H^H^nfWlj'Hn^ej (29) 
From Eq. (29) the tap coefficient G(n) is dependent on the value H^n) 
corresponding to the impulse response of the channel. In the case where 
5 temporal variations of the channel due to high transmission rate or the like, 
that is, where the duration of a received burst path is so short that the channel 
characteristic of the received signal does not change in that time interval, the 
tap coefficient G(n) is substantially constant. Accordingly, by precalculating 
the tap coefficient G(n) prior to the processing of the information symbol 
|* 10 period and using the precalculated value for processing of each information 

is! 

j J symbol, the computational complexity can be reduced as compared with that 

Jjj in the conventional system that derives the tap coefficient G(n) for each 

v3 information symbol. Moreover, since the above-mentioned temporal 

; variations need not be taken into account as well in the calculation of the tap 

I y 15 coefficient by Eq.(20), an adaptive algorithm can be used-this avoids the 

(0 necessity for conducting the inverse matrix calculation to derive the tap 

S3 

I a coefficient G(n), permitting reduction of computational complexity. 

In view of the above, in the first embodiment of the present invention, 
the tap coefficient calculation is conducted by an adaptive algorithm to obtain 
20 G(n) in a known signal period or training signal period prior to the processing 
of the information symbol period, and the value thus obtained is used to 
perform the processing in the information symbol period. This eliminates 
the needs for the inverse matrix calculation and the derivation of G(n) at each 
point in time of the information symbol, enabling the computational 
25 complexity to be made lower than in the conventional scheme. 

Figs. 3 and 4 depict examples of functional configurations for first 
processing and the second and subsequent processing for equalization by the 
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adaptive equalizer according to the first embodiment of the present invention. 
The illustrated adaptive equalizer corresponds to the equalizer 21 in Fig. 1. 
Fig, 3 shows the functional configuration for the first equalization processing, 
and Fig. 4 for the second and subsequent equalization processing. The 
5 broken line in Fig. 4 indicates the portion that is not used at the current point 
in time (in second and subsequent equalization processing) but used in the 
first equalization processing. 

The adaptive equalizer of the first embodiment comprises an impulse 
response estimating part 1 00, a linear filter part 1 10, a soft decision value 

10 calculating part 120, a switch 101, a likelihood calculating part 130, and a 
switch 160. The impulse response estimating part 100 is made up of a 
transversal filter 103, a subtracter 105, and an adaptive estimating part 106. 
The linear filter part 1 10 is composed of a transversal filter 1 11, a tap 
coefficient calculating part 112, and a subtracter 113. 

1 5 A description will be given, with reference Fig. 3, of the first 

equalization processing in the iterative equalization. In this instance, the 
switch 101 is connected to a known signal, for example, training signal input 
terminal 102, and the switch 160 is connected to an input terminal 104. In 
the impulse response estimating part 100, a training signal sequence b(n), 

20 provided as a known pattern from the training signal input terminal 102, and a 
difference value (vector) Rc(n), obtained by subtracting an output replica 
H m (n)b(n) of the transversal filter 103 by the subtracter 105 from a received 
signal sample value sequence R(n) from the input terminal 104 are used to 
estimate the impulse response H m (n) of the channel in a calculating part 106A 

25 of the adaptive estimating part 1 06. This impulse response estimating 

calculation needs only to minimize the square of Rc(n) through the use of an 
adaptive algorithm utilized in an ordinary adaptive equalizer. The estimated 
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impulse response value H m (n) having converged at the point of termination of 
the training signal sequence b(n) is stored in a storage part 106B. 

In the linear filter part 1 10, the received signal sample value sequence 
R(n) is input to the transversal filter 1 1 1 as a linear filter, by which it is 
5 subjected to linear filtering with the tap coefficient G(n) calculated in the tap 
coefficient calculating part 112, then the filtered output Z(n)=G H (n) H R(n) is 
subtracted by the subtracter 113 from the training signal sequence b(n), and 
the difference is input to the tap coefficient calculating part 1 12 to update the 
tap coefficient G(n) of the linear filter. This updating calculation utilizes an 
1 0 adaptive algorithm. That is, in the first iteration of equalization, since 
R(n)=Rc(n) from Eq. (27), the tap coefficient G{n) is so determined by the 
adaptive algorithm as to minimize ||b(n)-G(n) H R(n)|| from Eq. 16. On this 
account, the tap coefficient G(n) can be calculated with less complexity than 
in the case of conducting the inverse matrix calculation by Eq. (20). The 
1 5 updating of the tap coefficient G(n) stops with the termination of the training 
signal sequence b(n), and the tap coefficient at this time is held unchanged. 

The received signal after termination of the training signal sequence 
b(n), that is, the information-symbol received signal sample sequence R(n), is 
subjected to linear filtering by the transversal filter 1 1 1 with the most recently 
20 updated tap coefficient G(n). The filtering output Z(n) and the ultimately 
estimated channel impulse response value H m (n) are used to calculate, as the 
first equalization output, the soft decision value Xi[b(n)] in the soft decision 
value calculating part 120. In the first iteration of equalization, since the a 
priori information A, 2 [b(n)]=0, A 1 =[b(n)]=^ 1 [b(n)]. The soft decision value 
25 A,i[b(n)] is fed to the de-interleaver 22 (see Fig. 1). 

Next, a description will be given first of the principle of equalization 
in the second and subsequent rounds of iterative equalization and then of 
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concrete equalization processing with reference to Fig. 4. In the second and 
subsequent rounds of equalization processing, the received signal R(n) used 
in the first iteration of processing is input again. On the other hand, when 
the likelihood b'(n) is calculated by Eq. (12) from the output of the decoder 23 
5 (Fig. 1) produced for the result of the first iteration equalization, to be more 
precise, from the soft decision value A, 2 [b(n)] of the code bit fed from the 
interleaver 25, the likelihood b'(n) is approximately +1 or -1 if Eb/No (Bit 
energy-noise ratio) is high. In other words, the estimated intersymbol 
interference code vector of Eq, (13) can be approximated as expressed by the 
13 10 following equation. 

|2 B'(nM±l,±l, ±1,0, ±1,..., ±1, ±1] T (30) 

*,? 

?fj The use of thus approximated b'(n) permits approximation of Eq. (19) to the 

j " following equation. 

! , A(n) =Cov{B(n>BXn)}=diag{0,0,..., 0,1,0, ...,0,0} (31) 

?« 

[ y 15 Substituting Eq. (3 1 ) into Eq. (20), the tap coefficient G(n) of the linear filter 
|8 can be obtained by conducting the following calculation. 

S G(n)=[H L (n)H L H (n) WiJ^Hl (32) 

where H L (n) is the vector of one column of M by J rows expressed by 
H L (n)=[h(n; J-l)...h(n; 0)] T =[h 0 (n; J-l)...h M -i(n; J-l) 
20 ho(n; J-2). . .liM-i(n; J-2). . .ho(n; 0). . .h M -i(n; 0)] T . 

In this instance, the use of the Matrix Inversion Lemma in the 
calculation of the tap coefficient G(n) of the linear filter by Eq. (32) allows 
reduction of the computational complexity therefor. With such 
approximation of the estimated intersymbol interference code vector B(n) as 
25 expressed by Eq. (30), H L (n) becomes the vector of one column of M by J 
rows as mentioned above —this avoids the need for conducting the inverse 
matrix calculation, permitting substantial reduction of the computational 
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complexity as compared with that needed by the conventional scheme. Thus, 
according to the present invention, in the second and subsequent rounds of 
iterative equalization the computational complexity in the tap coefficient 
calculating part for deriving the tap coefficient G(n) of the linear filter can be 
5 reduced through utilization of the approximation by Eq. (3 1) and the Matrix 
Inversion Lemma. 

Fig. 4 illustrates in block form an example of the functional 
configuration in the second and subsequent rounds of equalization by the 
adaptive equalizer of the first embodiment. The change-over switch 101 is 
\t t 10 connected to the output of the likelihood calculating part 130, and the switch 
12 1 60 is connected to the output of the subtracter 1 05. The soft decision value 

ii ^2[b(n)] of each bit of the code sequence, which is the output from the 

jjj decoder 23 (the output from the interleaver 25) in the previous round of 
• iterative equation, is input to the likelihood calculating part 130 from a 

IU 15 terminal 130, which conducts the calculation of Eq. (12), outputting the 
W likelihood b'(n). The likelihood b'(n) is filter-processed by the transversal 
U filter 103 using the estimated channel impulse response value H m (n) stored in 

the storage part 106B, obtaining H m (n)B'(n) in Eq. (14) as a replica of the 
intersymbol interference in the received signal R(n). The output replica of 
20 the transversal filter 103 is subtracted by the subtracter 105 from the received 
signal value sequence R(n) to derive the difference value Rc(n) of Eq. (14). 

The tap coefficient G(n) of the linear filter is calculated by Eq. (32) 
using, as the input therefor, only the estimated channel channel impulse 
response value H L (n) stored in the storage part 106B unlike in the scheme for 
25 deriving the tap coefficient of the linear filter in the first iteration of iteration. 

This calculation utilizes the afore-mentioned Matrix Inversion Lemma. 
That is, when this Lemma is expressed by A=B~ , +CD~ 1 C H , where A and B 
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represent (M by M) positive definite matrixes, C an (M by N) matrix, D an (N 
by N) positive definite matrix, an inverse matrix of A is given by 

A-^B-BCCD+C^BBC)- 1 ^ (33) 
Applying this Lemma to the calculation of the inverse matrix in Eq. (32), 
5 Hl^H^^^B-'+CD-'c", 

Hl^HlV^D^C", 

o 2 I=B'\ 

h(n) =C, 

OH 

|j 10 H L H (n)=C H 

\i By using them to calculate Eq. (33), the inverse matrix calculation in Eq. (32) 

1^ is solved. Incidentally, Eq. (33) also contains an inverse matrix calculation 

u (D+C H BBC) _1 , but this inverse matrix is also scalar, and hence it can similarly 

! . be calculated. 

J *f 15 The thus calculated tap coefficient value G(n) is used to perform 

\% linear filtering of the difference value Rc(n) by the transversal filter 111, and 

I A the filtering result Z(n) and the estimated channel impulse response value 

H L (n) stored in the storage part 106B are used to calculate Eq. (25) in the soft 
decision value calculating part 120 to derive the soft decision value Ai[b(n)]. 
20 The equivalent amplitude \i(n) in Eq. (25) is given by the following equation 
by substituting Eq. (31) for A(n) in Fig. (23). 

M(n) = [H L (n) H [H L (n)H L (n) H +oVH(n)] 
-H L (n) H G(n) (34) 
Incidentally, in the first iteration of equalization the tap coefficient calculating 
25 part 112 of the linear filter part 110 may also use the Matrix Inversion 
Lemma to calculate the tap coefficient G(n). 

As described above, this embodiment permits reduction of the 
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computational complexity for deriving the tap coefficient of the linear 
equalizer that is used for iterative equalization. Fig. 5 shows an error rate 
characteristic obtained by computer simulations, for example, when an RLS 
(Recursive Least Squares) algorithm of excellent convergence characteristic is 
5 used as an adaptive algorithm for deriving the tap coefficient of a linear 
equalizer used for iterative equalization. 

In Fig. 5, the ordinate represents BER (Bit Error Rate) and the 
abscissa Eb/No (bit energy-noise ratio). The modulation system used was 
BPSK (two-phase shift), the information transmission rate was 12 Mbps, and 
□ 10 the maximum Doppler frequency indicating the fading variation rate was 
si 1 000 Hz. The frame configuration used is a unique word of a known pattern 

i q for first 128 symbols, then the channel impulse response H m (n) was estimated, 

1 2 and 128 symbols were added as information symbols to the end of the training 

j A symbols. The RLS algorithm was used for the estimation of the impulse 

f J 15 response. The channel was set as a five-wave Rayleigh fading channel of 

\i the same level, and the receiving side was to perform two-branch diversity 

13 

f * and iterative equalization by a five-tap linear filter. Further, SOVA 

(Soft-Output Viterbi Algorithm) was used for the decoding in the iterative 
equalization, and the number of iterations was two. 

20 In Fig. 5, the curves AO and Al show tap coefficients of the linear 

filter obtained by conducting the inverse matrix calculation at respective point 
in time in the cases where the number of iterations were one and two, 
respectively. The curves B0 and Bl show the tap coefficients G(n) derived 
by the scheme according to the first embodiment of the present invention, the 

25 curve B0 indicating the tap coefficient values by using the adaptive algorithm 
and the curve Bl the tap coefficient values by using the Matrix Inversion 
Lemma when the number of iterations was two. In the case of the present 
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invention, degradation of the bit error ratio is only on the order of 0.1 dB as 
compared with the ratio in the case of conducting the inverse matrix 
calculation at each point, it can be regarded as about the same as in the latter 
case. The present invention does not need to conduct the inverse matrix 
5 calculation expressed by Eq. (26), and hence it permits substantial reduction 
of the computational complexity for deriving the tap coefficient of the linear 
equalizer that is used for iterative equalization. Incidentally, it was 
confirmed that the same goes for the cases where the known pattern is 
128-symbol and the number of information symbols is 256 and where the 
13 10 known pattern is 1 28-symbol and the number of information symbols is 5 1 2. 
|I This means that when the duration of a burst path (the frame length) is short 

?fj about the same as these information symbol durations, impulse response 

jS variations of the channel can be ignored, and the bit error rate is 

L approximately equal to that in the case of conducting the inverse matrix 

J [j 15 calculation at each symbol time, and with less computational besides. 
SSJ Letting M and J represent the number of antennas (channels) and the length of 

intersymbol interference (the number of taps of the linear filter 111), 
respectively, the computational complexity is around 0((MJ) 3 ) in the case of 
conducting the inverse matrix calculation at each symbol time, but in the 
20 above-described embodiment of the present invention the computational 

complexity is approximately 0((MJ) 3 ), which is substantially lower than in the 
above case. 

While the examples of Figs. 3 and 4 are shown to include two 
transversal filters 103 and 1 1 1 with a view to facilitating a better 
25 understanding of the function configuration, the number of filters need not 
always be two but may also be one. Moreover, the tap coefficient 
calculating part 112 and the adaptive estimation part 106 are shown as 
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different blocks, but they can be combined into one unit 

As described previously, the one mode of the first embodiment utilizes 
the fact that the impulse response of the channel does not so vary over a short 
range of time, and hence it is characterized in that the training signal at the 
5 beginning of the received burst path (frame) or the known signal is used to 
obtain the estimated impulse response value H m (n), which is used to the 
subsequent processing. 

Further, according to the other mode of this embodiment, in the 
second iterations of iteration, only a J-th element of a diagonal matrix is set to 
10 1 which is a co variance matrix of the difference between the code bit vector 
B(n) and the estimated intersymbol interference code vector B'(n) as 
expressed by Eq. (3 1), and the other diagonal elements are approximated to 
zeros, by which the calculation for the tap coefficient G(n) becomes a vector 
calculation as indicated by Eq. (32)-- this avoids the necessity of conducting 
15 the inverse matrix calculation and hence allows reduction of computational 
complexity. In the first processing the known signal (training signal) need 
not always be used, but instead the output from a discriminating part 150, 
which discriminates the output from the transversal filter 103 between two 
value in binary, may be used in place of the signal from the output of the 
20 switch 1 01 as indicated by the broken lines in Fig. 3. The estimated impulse 
response value H m (n) ultimately obtained by the first processing is used in the 
second and subsequent processing. 

EMBODIMENT 2 

25 In the first embodiment described above with reference to Figs. 3 and 

4, all those of path components so-called multipath received signal R(n) 
arriving at the input terminal 104 via different paths which are within the 
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range of the maximum number J of delayed symbols to be considered are 
equalized as will be seen from Eqs. (5), (6), (8) to (11), (15) and (25). 
According to the channel environment of the channel, however, path 
components of relatively large received power may sometimes arrive apart in 
5 time from other path components as indicated by multipath components dl , 
d2, . in Fig. 6. In such a channel environment, too, the first embodiment 
performs equalization of all the received path components concerned. 

That is, in the first embodiment, path components within the range of 
the maximum number J of delayed symbols, even if very small, in extreme 
J 3 10 cases, zero in received power, are also subjected to equalization, and 
\2 consequently, much time is needed for calculation. As a solution to this 

| jj problem, this embodiment is adapted to omit particular ones of the 

fj dimensions corresponding to multipath components of negligibly small 

* received power to reduce the sizes of vectors and matrixes to be handled, 

|U 15 thereby reducing the computational complexity and hence shortening the 

§* 

id calculating time. 

ij 

\* Calculating the likelihood b'(n) of the soft decision value X 2 [b(n)] for 

the code bit string b(n) that is the decoder output for the result Ai[b(n)] of the 
first iteration of equalization, by Eq. (12) as in the second iterations of 

20 equalization in the first embodiment, the calculated value is virtually 4-1 or -1. 
Through utilization of this fact for the approximation of the estimated 
intersymbol interference code vector B'(n) as expressed by Eq. (30), the tap 
coefficient G(n) of the linear filter 11 1 can be calculated using the Matrix 
Inversion Lemma as follows: 

25 G(n)= [HL(n)H L (n) H + o^HkXn) 

- H L (n)/(a 2 + H L (n) H H L (n)) (35) 
where H L (n) is a column vector of a M by J length as shown 
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below. 

H L (n) » [ho(n; J-l)...h M -i(n; J-l) ho(n; J-2)...h M -i(n; J-2) ... 
...ho(n;0)...h M -i(n;0)] T (36) 
By the approximation of A(n) in Eq. (23) as indicated by Eq. (32), the 
5 equivalent amplitude u(n) of the output signal is also calculated as expressed 
by (34). 

On the other hand, Z(n) is given by 
Z(n) =G(n) H Rc(n) (37) 
as is the case with Eq. (15). By using in Eq. (25) the linear filter output Z(n) 
la 10 obtained by Eq. (37), the extrinsic information value Ai [b(n)] derived by the 
I* linear equalizer can be used as the results of second and subsequent 
iff equalizations. 

j 't Now, assume that path components j= 1 to j=J-2 within the range of the 

\x maximum number J of delayed symbols to be considered in the equalizer have 

SI 15 not arrived but that path components j=0 and j=J- 1 have arrived. If the 
SS impulse response of the channel can accurately be estimated, other path 

1 A components than effective ones are such that h(n; 1 ),h(n;2), . . . ,h(n; J-2)=0. 

Accordingly, it follows from Eq. (35) that 

H L (n)=[ho(n;J-l). . .hw-ifaJ-l) 0. . .0 0. . .0 ho(n;0) 

20 ...hM.,(n;0)] T (38) 

Since the tap coefficient G(n) of the linear filter can be calculated 
from H L (n) as shown by Eq. (35), it follows that 

G(n)= [go(n;J-l)...g M -i(n;M) 0...0 g 0 (n;0). . .g M -i(n;0)] T (39) 
The derivation of the soft decision value A^n)] by the linear 
25 equalizer by Eq. (25) requires the calculation of the equivalent amplitude u<n) 
of the output signal by Eq. (36) and the linear filter output Z(n) by Eq. (37). 
The equivalent amplitude ji.(n) and the linear filter output Z(n) are calculated 
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as shown below. 

^(n) = HL(n) H G(n) 

- ho(n; J-l)xgo(n;J-l)+...+ hM-i(n; J-l)xg M .,(n;J-l)+...+OxO+... 
. . .+0x0+. . .+ ho(n;0)xgo(n;0)+. . .+ h M -i(n;0)xg M . I (n;0) 

(40) 

Here, in the case where zero-elements of HiXn) H and G(n) are omitted, 
that is, where the dimensions of H L (n) H and G(n) are lowered as shown below, 
too, the value of u(n) naturally remains unchanged. 

H L '(n) H = tho(n;J-l)...h M -,(n;J-l) ho(n;0)...h M .,(n;0)] (41) 
G'(n) = [go(n;J-l)...g M -i(n;J-l) go(n;0)...g M -i(n;0)] T (42) 
u(n) = H L '(n) H G'(n) 

= ho(n;J-l)xgo(n;J-l)+...+ h M -i(n;J-l)xg M . 1 (n;J-l) 

+h 0 (n;0)xg 0 (n;0)+. . .+h M -i(n;0)xg M .,(n;0) (43) 
And, Z(n) becomes as follows: 
Z(n) = G(n) H R c (n) 

= go(n; J-l)xr 0 (n+J-l)+...+ g M .i(n; J-l)xr M .,(n+J-l)+... 

+0xr 0 (n+J-2)+. . .+0xr M .,(n+J-2)+ +0xr 0 (n+l)+. . .+0x 

r M -i(n+l)+go(n;0)xr 0 (n)+. . .+gM-i(n;0)xr M -i(n) (44) 
In Eq. (44), since G(n) H has O-elements, the values of Rc(n) 
corresponding to the O-elements are unnecessary in the derivation of Z(n) as 
shown below. 

^•(n) = [r 0 (n+J-l). . .r M -i(n+J-l)r 0 (n). . .r M -i(n)] (45) 
Z'(n) = G'(n) H R <: , (n) 

= go(n;J-l)x ro (n+J-l)+. . .+ gM-i(n; J-l)xr M .i(n+J-l)+. . . 
. . .+ g 0 (n;0)xr o (n)+. . .+ gM-i(n;0)xr M .,(n) (46) 
Since the dimension of Rc(n) can be lowered as described above, the 
dimension of the received signal vector R(n) can also be lowered by Eq. (14). 
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Further, since the dimension of the intersymbol interference component 
replica H ro (n)B'(n) can also be reduced, the dimension of the channel matrix 
H m (n) can be reduced as well. 

From the above, it can be seen that calculations for the path 
components having arrived can be dispensed with in the calculation for 
deriving the extrinsic information by the linear equalizer. That is, in the 
estimated impulse response value of the channel, the tap coefficient G(n) of 
the linear filter, the received signal vector R(n) and Rc(n) obtained by 
subtracting the intersymbol interference component from the received signal 
vector, those elements of vectors and matrixes corresponding to path 
components having not arrived can be removed. Accordingly, it is possible 
to decrease the number of dimensions of vector components and the number 
of elements of matrixes necessary for the calculation involved in the 
derivation of the soft decision value of the linear equalizer. 

However, the estimated impulse response value does not become zero 
due to an error in the estimation of the impulse response of the channel. 
Hence, a certain threshold value ha, is set for comparison and discrimination 
of the estimated impulse response value as described below. For example, if 
h(n;l),h(n;2),...h(n;J-2)<h th , then it is regarded as 
h(n; 1 ),h(n;2), . . . h(n; J-2)=hth. Letting L represent the number of effective 
received path components larger than the threshold value hth, the dimension of 
the estimated impulse response vector H L (n), which indicates all estimated 
impulse response values within the range of the maximum number J of 
delayed symbols, can be lowered from an M*J-th dimensional vector to an 
M*L-th dimensional vector H' L (n) as expressed below. 

Hi/n) =[h 0 (n; J-l)...h M -i(n; J-l) 0 ... 0 ho(n;0)...hM-i(n;0)] T (47) 
H L '(n)=[ho(n; J-l)...h M -i(n; J-l)h o (n;0)...h M .,(n;0)] T (48) 
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Further, the number of row elements of the channel matrix H„,(n) can 
be lowered from 2J-1 to 2L-1. 



H m (n) = 



h(n;0) h(n;l) - h(n;J-l) 0 
0 h(n;0) h(n;l) •• h(n;J-l) 0 



0 
0 



H, 



m 



0 h(n;0) h(n;l) 
0 h(n;0) 



••• h(n;J-l) 0 
0 b(n;0) 0 



h(n;J-l) 0 
h(n;l) ••• h(n;J-l) 



(49) 



0 1 

h(n;J-l)J 



MLx(2J-l) 



'JMJx(2J-l) 



(50) 



If 



h(n; l),h(n;2),. . .h(n;J-2)^w 



then 



go(n; 1). . .gM-i(n; l),go(n;2). . .g M -i(n;2),. . .,go(n;J-2). . .g M -i(n;J-2)s 0 
1 0 is set; the tap coefficient vector G(n) of the linear filter in the second 
iterations of iteration is given by the following equation. 

G'(n)= [go(n;J-l)...gM-,(n;J-l) go(n;0)...g M .i(n;0)] T (51) 
This is apparent from Eq. (3 5) in the calculation of H L '(n) H H L '(n). Moreover, 
the tap coefficient vector G(n) can also be lowered from the M* J-th 
1 5 dimensional vector to the M*L-th dimensional vector G'(n). 

Similarly, the received signal vector R(n) and the difference vector 
R^n) can also be lowered from the M* J-th to the M*L-th dimension. 

The equivalent amplitude u(n) of the output signal can be calculated 
by Eq. (43) from H L '(n) and G'(n), and the linear filter output Z'(n) can be 
20 calculated by Eq. (46) from G'(n) H and Rc'(n). Accordingly, they can each 
be derived by calculations of vectors of low dimensionality. 

Furthermore, it is possible to lower the dimension of the number of 
column elements at positions of path components regarded as zeros in the 



channel matrix Hm'(n) of Eq. (50) generated from the estimated impulse 
response value as expressed by the following equation, and the dimension of 
the channel matrix for deriving the soft decision value by Eq. (25) can also be 
reduced. In an example in which only paths of j=0 and j=J are set as 
effective paths, the number of column elements can be decreased from (2J-1) 
to (2L-1). 

t»o> 0 J, J < 52 > 

L 0 h ( D >°) h ( n ' J " l )JMLx(2L4) 

Figs, 7 and 8 illustrate in block form examples of functional 
configurations in the first iteration and second and subsequent rounds of 
equalization. In this embodiment, an identifying part 140 and a switch 170 
are additionally provided in the functional configuration of the first 
embodiment shown in Figs. 3 and 4, and the tap coefficient calculating part 
1 1 2 is composed of a calculating part 1 1 2A, a storage part 1 1 2B and an 
identifying part 1 12C. The adaptive estimation part 106 is also added with 
an identifying part 106C. 

In the first place, the switch 101 is connected to the terminal 102, and 
the switches 160 and 170 are both connected to the terminal 104. In the first 
iteration of equalization the processes for obtaining the estimated impulse 
response value H m (n) and the tap coefficient G(n) are the same as in the first 
embodiment. Further, in Fig. 7, effective paths are identified by the 
identifying part 106C based on the estimated impulse response value H m (n) 
having converged at the end of the training signal. Assume, for example, 
that estimated impulse response values are obtained at respective symbol 
delay times within the range of the maximum number J of delayed symbols as 
shown in Fig. 9. The effective path identifying method is, for instance, to 
select paths above a certain threshold value h^ as depicted in Fig. 10. The 
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threshold value hth is for instance, the average power of all received paths, and 
received paths above the average power are selected The effective paths 
thus identified are each represented, for example, by the number of delayed 
symbols and stored in the storage part 106B together with the corresponding 
5 tap coefficient value h. This data may be stored as data of the vector H L '(N) 
given by Eq. (48) or as data of the matrix H m '(n) given by Eq. (48); and it can 
easily be obtained from the one to the other and vice versa. Further, only 
those of the tap coefficient vectors G(n) of the linear filter 1 1 1 which 
correspond to the effective paths identified by the identifying part 106C are 
C3 1 0 each stored as the tap coefficient vector G'(n) in the storage part 1 12B. 
I~ Moreover, the switch 160 is switched to the identifying part 140, in 

IH which only an information symbol received signal R'(n) corresponding to the 

I* effective path in the signals received after the training signal, that is, in the 

§ * information symbol received signals R(n), is selected, and the thus selected 

IU 

U 15 signal is subjected to linear filtering by the transversal filter 111 using its tap 

iS 

ig coefficient G f (n) stored in the storage part 1 12B most recently. The linear 

filtering output Z'(n) and the most recently estimated impulse response value 
H' L (n) are used to calculate the soft decision value Aj[b(n)] in the soft 
decision value calculating part 120. 

20 In the second and subsequent rounds of equalization, the switch 101 is 

connected to the likelihood calculating part 130, the switch 160 to the 
identifying part 140 and the switch 170 to the output of the subtracter 105 as 
depicted in Fig. 8. The likelihood b'(n) representing the estimated 
intersymbol interference code, provided from the likelihood calculating part 

25 130, is subjected to filtering by the transversal filter 103 using the estimated 
impulse response value H L (n) stored in the storage part 106B, and the output 
replica H L (n)b'(n) is subtracted by the subtracter 105 from the received signal 
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sample value sequence R'(a) to obtain the difference value Rc'( n )- For only 
the path identified in the identifying part 106C of the impulse response 
estimating part 100, the difference value R</( n ) is identified by the identifying 
part 140, and the identification result is input to the linear filter 111. The tap 
5 coefficient G'(n) of the linear filter 1 1 1 is derived from only the estimated 
impulse response value H L '(n) stored in the storage part 106B, by Eq. (35) 
through utilization of the Matrix Inversion Lemma, and the tap coefficient 
G'(n) is stored in the storage part 1 12B. The thus stored estimated impulse 
response value H L '(a)> the tap coefficient G'(n) of the linear filter 111 and the 

13 10 identified difference value R</( n ) are used to calculate the soft decision value 

(3 

I* Ai[b(n)] in the soft decision value calculating part 120. 
in As described above, the second embodiment permits reduction of the 

\& computational complexity for deriving the soft decision value of the linear 

§* equalizer. Fig. 1 1 shows the results of computer simulations of the bit error 

u 15 rate characteristic in the case of using the RLS (Recursive Least Squares) 
q algorithm of excellent convergence characteristic for deriving the tap 

coefficient of the linear equalizer in the first iteration of iterative equalization 
and the Matrix Inversion Lemma in the second and subsequent iterations of 
iterative equalization. 
20 The modulation system used is BPSK; the information transmission 

rate is 10 Mbps; the maximum Doppler frequency indicating the fading 
variation is 1000 Hz; and the frame configuration is one that a training signal 
of a known pattern of 63 symbols is added as an information signal to 256 
symbols. A channel, in which only path components j=0 and j=4 as 
25 indicated by the solid lines in Fig. 12 had arrived, was assumed as a channel 
model. The receiving side performs two-branch diversity and uses SOVA 
(Soft-Output Viterbi Algorithm) for decoding. On the assumption that the 
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channel was ideally estimated, the number of iterations was set to three. 

In Fig. 1 1 , the curves Al, A2 and A3 show the results obtained in the 
cases where: equalization was performed once and iterated twice and three 
times, respectively; the tap number of the linear equalizer was considered for 
5 0 to 4T delayed paths; and the entire information of five taps was used as the 
tap coefficient of the linear equalizer. The curves Bl , B2 and B3 show the 
results obtained in the cases where: equalization was similarly performed 
once and iterated twice and three times, respectively; and two-tap (only 0 and 
4T delayed paths) information was used as the tap number of the linear 
10 equalizer. 

As is evident from Fig. 11, a comparison of the results obtained with 
the use of five-tap and two-tap information as the tap number of the linear 
equalizer indicates that two and three iterations of equalization provide about 
the same error rate characteristic with low computational complexity. 

15 

EFFECT OF THE INVENTION 

As described above, the first embodiment of the present invention 
enables the linear equalizer for iterative equalization enables to reduce the 
computational complexity for calculating the tap coefficient of the linear filter 

20 and hence shorten the time for the calculation. 

According to the second embodiment of the invention, in a channel 
environment in which path components of relatively large power have arrived 
apart in time from other path components, it is possible to adaptively reduce 
the computational complexity for calculating the extrinsic information by the 

25 linear equalizer that performs iterative equalization. 
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